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Different phenomenological RG transformations based on scaling relations for the derivatives of the inverse 
correlation length and singular part of the free-energy density are considered. These transformations are tested 
on the 2D square Ising and Potts models as well as on the 3D simple-cubic Ising model. Variants of RG equations 
yielding more accurate results than Nightingale's RG scheme are obtained. In the 2D case the finite-size equations 
which give the exact values of the critical point or the critical exponent are found. 



1. Introduction 

The phenomenological renormalization-group 
(RG) method [jjj is a powerful tool for the in- 
vestigation of critical phenomena. As it is known 
H), phenomenological RG can be constructed by 
using not just the correlation length as it is done 
in Nigtingale's approach but using any other 
quantity with a power-law divergence at criti- 
cality. Binder || suggested a phenomenologi- 
cal renormalization scheme by using the order- 
parameter moments (cumulants) which are, on 
the one hand, related to the higher susceptibilities 
and, on the other hand, immediately suitable for 
the Monte Carlo simulations. Recently Itakura 
jd] extended Binder's cumulant crossing method 
taking linear combination of several different re- 
duced moments. 

In this report, I discuss various RG transfor- 
mations which follow from general scaling func- 
tional equations. These equations are evaluated 
in terms of the eigenvalues and eigenvectors of 
the transfer matrices. By large transverse sizes 
of partly finite subsystems, all those transfor- 
mations must yield the same results. However, 
for the small sizes which normally are used in 
practice, different RG equations lead to estimates 
with distinct accuracies. My aim is to find the 



best strategies of a phenomenological renormal- 
ization group method. This is especially impor- 
tant for 3D systems. 

2. Phenomenological RG Equations 

Let us write the finite-size scaling equations for 
the derivatives of the inverse correlation length 
kl and the singular part of the reduced free- 
energy density /£: 
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where z^ m,n \x,y) = d m+n z/dx m dy n (z is kl or 
/£), t = K — K c is the deviation from critical cou- 
pling, h is a normalized external field, yt and yh 
are, respectively, thermal and magnetic critical 
exponents of the system, d is the space dimen- 
sionality, L is a characteristic size of a subsystem 
and b = L/L' is the rescaling factor. 

In the phenomenological approach proposed by 
Nightingale Q], eq.(0) with m = n = is com- 
bined with the ordinary expression for the inverse 
correlation length 

K L = H\[ L) /xi L) ), 



(3) 



in which \[ L ^ and X^' are the largest and second- 
largest eigenvalues, respectively, of the associated 
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transfer matrix. In the absence of a symmetry 
breaking field, the critical coupling K c is esti- 
mated from the equation 

Lk l (K c ) = {L-1)k l _ 1 (K c ). (4) 

In writing this equation, one sets L' = L — 1. 

Another possible way to produce a phenomeno- 
logical renormalization group is obtained by using 
eq.(||) with m = n = 0. The fixed point is given 
by the relation 

L d fl(K c ) = (L-l) d fl_ 1 (K c ). (5) 

The dimensionless free-energy density, /l = foo + 
/£, of a subsystem L^ 1 x oo is calculated by the 
formula 



Table 1 



f L = L l - d \n\[ 



(L) 



(G) 



and the "background" foo is introduced as an ex- 
tra parameter. 

Besides eqs.(||) and (||), in this paper I also con- 
sider the following RG equations (resulting from 
the relations (|lj) and (||)): 



(4) 
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where X l = & 2 h/ dK 2 \ h=Q = f L {0 > 2} (K,0) 
is the zero-field susceptibility and Xl = 

| 9 4 /i/5^ 4 U=o = fi} 0,i \^j 0) 1S a nonlinear sus- 
ceptibility (eq.(|^) corresponds to Binder's phe- 
nomenological renormalization group); 

i 2 - d (4 1) ) 2 /xL = (i-i) 2 
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/XL = {L-l) 



l/X 



L-l, 



L^ d Kf/xl = {L-l) 1 - 2d ^l 1 /x 



KV >{K,0). 
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(10) 
Ex- 



Here, 4 n) = d n K L /dh n \ h=0 = 4 0,n) 

pressions for the derivatives of the inverse corre- 
lation length and the free energy with respect to 
h in terms of eigenvalues and eigenvectors of the 
transfer matrix are available in p[. 

3. Results and discussions 

To represent numerical data in tables, eqs.([|), 
§)> ©> (!) and © will be labeled by sym- 



and 
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respectively. 



Estimates 


of iiTc for 


the 2D sq Ising 


lattice; 


j^exact _ 0.440 686 . . . 








eq. 


(2,3) 


(3,4) 




(4,5) 


(«) 


0.42236 


0.43088 (- 


-2.23%) 


0.43595 


(x (4) /x 2 ) 
(« (4) /x 2 ) 


0.42593 


0.43242 (- 


-1.88%) 


0.43672 


0.42596 


0.43243 (- 


-1.87%) 


0.43673 


(/ s ) 


0.44324 


0.44168 (- 


-0.23%) 


0.44105 


( K ( 2 )/ x ) 


0.47420 


0.45153 (- 


-2.64%) 


0.44626 



|ln2 



In table [l], results for the critical coupling in the 
Ising model on a square lattice are given. The cal- 
culations were carried out for strips L x oo with 
a periodic boundary condition in the transverse 
direction. The estimates are shown for the pairs 
(L — 1, L) with L < 5. In the case of (3, 4) pairs, 
the errors are also given. The type of phenomeno- 
logical RG equations which have been used are 
indicated in the first column of the table. In this 
model K c = \ ln(l + V2) and = 2G/vr 
(G is Catalan's constant) p. 

It is seen from table fythat the best lower 
bound is given by eq.([10|). Slightly worse re- 
sults are obtained by Binder's phenomenological 
renormalization-group procedure. This approach 
which is normally implemented by Monte Carlo 
simulations was used in the transfer-matrix ver- 
sion in Nightingale's renormalization (first 
line in table 1) which is traditionally used by 
transfer-matrix calculations has only the third 
position in accuracy among the lower estimates. 

I also found the phenomenological RG equa- 
tions leading to the upper bounds for K c (last 
two lines in table 0). Among these more accurate 
results are provided by eq.(^). The magnitude 
of the error in line four of table 1 is the least 
among all lower and upper estimates of K c . Un- 
fortunately, such approach requires a knowledge 
about the background f^. 

Transfer-matrix eigenvalues for the Ising strips 
are known in analytical form ||. Using this fact 
I considered the RG transformation (||) which 
makes use of the first derivative with respect to 
K (m = 1, n = 0; h = 0). For the fixed point this 
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Table 2 

Estimates of K c for the 2D sq 3-state Potts lat- 
tice; Kf xact = 1.005 052... 



cq. 



(2,3) (3,4) 



(4,5) 



(k) 0.96248 0.98350 (-2.1%) 0.99467 
(/t (1) 7x) 0.99311 0.99920 (-0.6%) 1.00380 
(f s ) 1.00927 1.00667 (+0.2%) 1.00565 



transformation gives 

L d -y*u s L (K c ) = (L- if-ytul^iKc), (11) 

where u s L — ul — is the singular part of the 
reduced energy density, and ul — dfhjdK. Re- 
markably, the root of eq.(O) is equal to the ex- 



act value of K c since u s L {K c ) = for all L. In 
other words, all finite-size corrections to the back- 
ground u oc (= %/2 Q) are zero. 

Moreover, in the 2D Ising model dn^/dK at 
K = K c is also independed of L and, therefore, 
eq. ([!]) with m = 1 and n = gives the exact value 
for the critical exponent: v = 1/yt = 1- 

The data obtained for the 3-state square Potts 
lattice are collected in table |^. For this model 
K c = ln(l + a/3) and /«, = 4G/3tt + bi(2v / 3) + 
| ln(2 + \/3) J|. Inspecting table ||, it is seen that 
eqs.(||) and (||) lead to more qualitative estimates 
than Nightingale's approach. Again, the lowest 
absolute error is yielded by the phenomenological 
RG equation based on /£. 

Numerical calculations on strips L x oo show 
us that in the 2D g-state Potts model the finite- 
size corrections to the background energy, Uoo = 
1 + lf-y/q are also absent and, consequently, 
the equation 



u L (K c ) = u L '(K c ) 



(12) 



yields the exact value of K c . Note that this equa- 
tion has been derived earlier from other consider- 
ations ||. 

Let us discuss now the results presented in ta- 
ble H for the 3D Ising model on a simple-cubic lat- 
tice. For this model K c = 0.221 6544(3) ffij and 
/oo = 0.77790(2) jn]]. Renormalizations were 
done for the L x L x oo parallelepipeds with pe- 
riodic boundary conditions in both transverse di- 



Table 3 

Estimates of K c for the 3D sc Ising lattice; 
K exact _ 0.221 6544(6) 



eq. 



(2,3) (3,4) 



(x (4) /x 2 ) 
(* (4) /x 2 ) 

(/ s ) 



0.21340 0.21826 (-1.53%) 

0.21823 0.22002 (-0.74%) 

0.21824 0.22006 (-0.72%) 

0.22354 0.22236 (+0.32%) 

0.22658 0.22314 (+0.67%) 



rections. As in the 2D case, the best lower values 
of K c are obtained from eq. (|To|) ■ 

In the 3D case the amplitudes of the finite-size 
corrections to the critical-point energy are not 
equal to zero. As a result, eq.(12) only yields 
an approximate value of K c . 
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